Strength of the d x 2_ y 2 pairing in the two-leg Hubbard ladder 
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In the ground state of the doped two-leg Hubbard ladder there are power-law decaying d x 2_ y 2- 
type pairing correlations. It is important to know the strength and the temperature scale of these 
correlations. For this purpose, we have performed determinantal Quantum Monte Carlo (QMC) cal- 
culations of the reducible particle-particle interaction in the Hubbard ladder. In this paper, we report 
on these calculations and show that, at sufficiently low temperatures, resonant particle-particle scat- 
tering takes place in the d x 2_ y 2 pairing channel for certain values of the model parameters. The 
QMC data presented here indicate that the d x 2_ y 2 pairing correlations are strong in the Hubbard 
ladder. 

PACS numbers: 71.10.Fd, 71.10.Li, 74.20.Rp 
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It is generally accepted that the high-T c cuprates Q 
are ^2^2 -wave superconductors It is also known 

that the spin-fluctuation exchange in electronic models 
0, |]| and an electron-phonon interaction which has its 
largest coupling at small wavevectors [E 01 can lead 
to an effective electron-electron interaction which is at- 
tractive in the d-w&ve Bardeen-Cooper-Schrieffer (BCS) 
pairing channel. However, it is not yet clear what is 
the pairing interaction responsible for superconductivity 
in the high-T c cuprates. Within this context, it is im- 
portant to determine the maximum possible strength of 
the d x 2_ y 2 pairing correlations which one can obtain in 
a purely electronic model such as the Hubbard model. 
Unfortunately, the ground state of the two-dimensional 
(2D) Hubbard model remains beyond the reach of the 
exact many-body techniques However, in the case of 
the two-leg Hubbard ladder, the Density Matrix Renor- 
malization Group (DMRG) calculations found power-law 
decaying d x 2 _ y 2 -type pair-field correlations 

[IB, EH- This 

is probably the only model where it is known from exact 
calculations that the pairing correlations get enhanced by 
turning on an onsite Coulomb repulsion in the ground 
state |ll|,ll2j. Hence, it is important to determine the 
maximum possible strength and the temperature scale of 
the d x 2_ y 2 pairing correlations in the Hubbard ladder, 
which is the main motivation of this paper. 

Here, we present Quantum Monte Carlo (QMC) results 
on the reducible particle-particle interaction T in the BCS 
channel, which is illustrated in Fig. 1. The reducible T 
serves as a powerful probe of the pairing correlations. For 
example, in the case of an s-wave superconductor, T at 
the Fermi surface would diverge to —00 due to repeated 
particle-particle scatterings in the BCS channel, when 
the superconducting transition is approached, T — > T c + . 
In this paper, we investigate the strength of the d x 2_ y 2 
pairing in the Hubbard ladder by making use of the exact 
QMC data on T. 

In the following, we show that the reducible particle- 
particle interaction T in the Hubbard ladder exhibits di- 
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FIG. 1: Feynman diagram for the reducible particle-particle 
interaction F(p'\p) in the BCS channel. Here, p denotes 
(p,iu) n ) with Matsubara frequency uo n . In this diagram, the 
incoming fermions at states p with up spin and — p with down 
spin scatter to states p' with up spin and —p' with down spin 
by exchanging q = p' — p. 



verging behavior as the temperature is lowered for cer- 
tain values of the model parameters. We find that, on the 
Fermi surface, T can become strongly repulsive (attrac- 
tive) for q « (n, it) (q ks 0) momentum transfers, which 
correspond to backward (forward) scatterings. In partic- 
ular, near half-filling and for Coulomb repulsion U = it 
and temperature T = O.lt, where t is the hopping ma- 
trix element, the backward and forward scattering am- 
plitudes can become an order of magnitude larger than 
the bare Coulomb repulsion or the bare bandwidth. This 
type of momentum dependence of T implies that reso- 
nant particle-particle scattering is taking place in the 
BCS channel already at T « O.lt. For U — 8t, we observe 
similar behavior at about twice higher temperatures. The 
temperatures studied in this paper are lower than those 
reached in previous QMC calculations of V for the 2D 
[l3| and the two-leg |14| Hubbard models. In addition, 
we present results on the solution of the Bethe-Salpeter 
equation in the BCS channel, which quantitatively de- 
termines the strength of the pairing correlations. The 
QMC data shown in this paper imply that the d x 2_ y 2- 
type pairing in the two-leg Hubbard ladder is strong for 
certain values of the model parameters. 

We begin by briefly describing the previous studies of 
d x 2_ y 2 pairing in the Hubbard ladder. The DMRG cal- 
culations found that the rung-rung pair-field correlation 
function decays as power law in the ground state of the 
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doped Hubbard ladder The mean-field calculations 
suggested the d^2_ ,.2 t ype of symmetry for pairing in 
doped spin ladders [lfij . The exact-diagonalization |l2j 
and the DMRG ^lj calculations found that the d x 2_ y 2 
pairing correlations are most enhanced when the inter- 
chain hopping *j_ is greater than the intrachain hopping 
t, in particular, for t± ks 1.5* in the intermediate cou- 
pling regime and near half-filling. In this case, the QMC 
calculations showed that the irreducible particle-particle 
interaction peaks for momentum transfers near (jr, tt) due 
to antiferromagnetic fluctuations [T^j . 

The two-leg Hubbard model is defined by 

i,\.<j i.a 

+h.c.) + ^i,A,T^i,A,| - A* X! n i,K°> C 1 ) 

i,A i,A,cr 

where * is the hopping parameter parallel to the chains 
(along x), and tj_ is for hopping perpendicular to the 
chains (along y). The operator c] A a (ci.\ t a) creates (an- 
nihilates) an electron of spin a at site i of chain A, and 
n>i,\,<7 = c| Ao .Ci.A, CT is the electron occupation number. 
As usual, U is the onsite Coulomb repulsion, and /i is 
the chemical potential. In addition, periodic boundary 
conditions were used along the chains, and t± denotes 
interchain hopping for open boundary conditions. 

In obtaining the data presented here, the determinan- 
tal QMC technique described in Ref. [16| was used. The 
calculation of the reducible interaction T follows the pro- 
cedure described in Ref. 0] for the 2D case. The BCS 
component of the reducible particle-particle interaction 
T(p'\p) is illustrated in Fig. 1, where p — (p,ico n ) with 
Matsubara frequency ui n = (2n+ \)ttT. In the following, 
results will be shown for the reducible interaction in the 
singlet channel, T s (p'\p) = \ [r(p'|p) p'\pj\ ■ In par- 

ticular, the momentum dependence of r s (p', iuj n t |p, iu) n ) 
will be shown at the lowest Matsubara frequency lu u — 
u> n i = ttT. In the ground state and for U — 4*, the 
d a .2_ y 2-type pairing correlations are most enhanced near 
half-filling for t± « 1.6* [ll|. When U = St, this oc- 
curs for t± w 1.4*. In this paper, the QMC data will 
be presented using these two parameter sets for a 2 x 16 
lattice. 

The momentum structure in T s (p'\p) depends sensi- 
tively on where p' and p are located with respect to the 
Fermi surface. Hence, we will first show results on the 
single-particle spectral weight A(p,u>) = — ^lmG(p,u> + 
iS), where G is the single-particle Green's function, and 
discuss the location of the Fermi-surface crossing points. 
We have obtained A(p, u) from the QMC data on the 
single-particle Green's function along the Matsubara- 
time axis by using the maximum-entropy analytic con- 
tinuation method. Figure 2 shows A(p, ui) versus u> for 
T = 0.1*, U = At, t x = 1.6* and (n) = 0.94. Here, 
we see that the Fermi level crossing for the antibond- 
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FIG. 2: (color online) Single-particle spectral weight A(p,u) 
versus ui at various p. The solid and dotted curves represent 
the bonding (p y = 0) and antibonding (p y = 7r) bands, re- 
spectively. These results are for T — O.lt, U = it, t± — 1.6* 
and (n) = 0.94. The arrows denote the quasiparticle positions 
for the U = case. 

ing (p y = 7r) band occurs for p x between ir/8 and 7r/4, 
while for the bonding (p y = 0) band there is spectral 
weight pinned near the Fermi level for 37r/4 < p x < it 
at this temperature. These results for Mp, w) are sim- 
ilar to those presented in Refs. [ill Il4] . In the fol- 
lowing, r s (p', iirT\p, mT) versus p will be shown for 
p' = (tt/4, 7t) near the Fermi level and for p' = (0,7r) 
at the saddle point. 

Figure 3(a) shows r s (p', inT\p, inT) versus p while 
p' is kept fixed at (7r/4,7r). In the left panel, p x is 
scanned from — 7r to 7r for p y = 0, while in the right 
panel p x is scanned for p y = n. Here, it is seen that 
repulsive and attractive peaks develop in T s , as T de- 
creases from 0.25* to 0.1*. In particular, in the left panel 
it is seen that, when p « (— 37r/4, 0), a repulsive peak 
develops in T s , which corresponds to a scattering pro- 
cess with momentum transfer q « (tt,tt). In addition, 
in the right panel it is observed that a dip develops in 
T s when p k p' — (tt/4, tt) corresponding to zero mo- 
mentum transfer. This dip is due to resonant scattering 
in the d 2 ,2_ y 2-wave BCS channel. In a three-dimensional 
infinite system, when a c? a; 2_j / 2-wave superconducting in- 
stability is approached, r s (p', «7rT|p, iirT) at the Fermi 
level diverges to +oo for backward scattering, and to -co 
for forward scattering, which will be further discussed be- 
low. In Fig. 3(a), it is seen that r s is developing this type 
of repulsive and attractive peaks at T of order 0.1*. 

Figure 3(b) shows r s (p', iirT\p, iirT) versus p while p' 
is kept fixed at the saddle point (0,7r). In this case, T s 
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FIG. 3: Reducible particle-particle interaction in the singlet 
channel r s (p', iirT\p, inT) versus p x . Here, p x is scanned for 
p y = (left panel) and for p y — tv (right panel). In (a), p' is 
kept fixed at (tt/4, 7r), while in (b), p' = (0, 7r). Here, T s is 
shown in units of t for U = 4t, fx = l-6t and (n) = 0.94. 



develops a peak when p = (±7r, 0), corresponding to scat- 
tering with q = (tt,tt) momentum transfer. The magni- 
tude of this peak is comparable to that seen in the left 
panel of Fig. 3(a). However, the behavior for q = 
momentum transfer is different. As observed in the right 
panel of Fig. 3(b), r s for q = scattering remains pinned 
near zero for T down to 0.125i, and becomes attractive 
only below this temperature. Hence, at the saddle point 
(0, tt), the resonant scattering in T s for q = momentum 
transfer is weaker compared to that at (7r/4,7r). 

Figure 4 shows the T dependence of the backward and 
forward scattering components of T s . Here, r s is plotted 
as a function of the inverse temperature (3 for momentum 
transfers q = (tt, 7t) and (tt/8, 0), while p' is kept fixed 
at (tt/4, tt). This figure shows that, near the Fermi level, 
the backward and forward scattering amplitudes increase 
rapidly at low T, becoming an order of magnitude larger 
than the bare Coulomb repulsion at T ss O.lt. 

Figures 3 and 4 display the main features of T s , which 
can be summarized as follows: At low T, there are 
strong q » (tt, tt) scatterings over the whole Brillouin 
zone, while the q f» scatterings are most attractive 
near the Fermi surface. These features were observed 
at (n) = 0.94 and 0.875, and for U = 4f and 8t. The 




FIG. 4: Reducible particle-particle interaction in the singlet 
channel r s (p', i7rT|p, inT) versus the inverse temperature /3 
for backward (q = (n,n)) and forward (q = (tt/8, 0)) scatter- 
ings. Here, p = p' + q and p' is kept fixed at (7r/4, tt). These 
results are for U = it, t± — 1.6t and (n) = 0.94. 



one-loop Renormalization-Group (RG) technique is also 
utilized to obtain the momentum dependence of T in the 
Hubbard model [rH fl8j |. It would be useful to make 
comparisons of the QMC and the RG results for the mo- 
mentum dependence of T in the ladder case. 

In order to determine the strength of the pairing cor- 
relations, the Bcthc-Salpctcr equation for the reducible 
particle-particle interaction in the BCS channel, 
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was solved for X a and the corresponding eigcnfunctions 
<fra(p,iu) n ) by using QMC data on T and the single- 
particle Green's function G(p). For a three-dimensional 
infinite system, when the maximum A Q reaches 1, this 
signals a BCS instability to a state where the pair wave 
function has the form of the corresponding eigcnfunc- 
tion. For a one-dimensional system, A Q 's will always be 
less than 1, however, the T dependence of A Q 's gives in- 
formation about the characteristic temperature scale of 
the pairing correlations. For instance, if the maximum 
eigenvalue reaches 0.9 at some temperature, then this 
means that the leading pairing correlations are enhanced 
by a factor of 10 through repeated particle-particle scat- 
terings in the BCS channel. Hence, at this temperature, 
the system would exhibit strong pairing fluctuations. 

At low temperatures, the maximum A a of the 
Bethe-Salpeter equation corresponds to an eigenfunction 
c/)d(p 1 iuj n ) which has d x 2_ y 2 type of symmetry in the 
sense that it changes sign as p goes from (tt, 0) to (0, tt) in 
the Brillouin zone [l4| . The temperature evolution of the 
d x 2_ y 2-w&ve irreducible eigenvalue is shown in Figure 
5(a) for U = 4t and t± = 1.6i. As T decreases, grows 
monotonically reaching 0.75 at T = O.lt for (n) = 0.94. 
The simple extrapolation of these results suggests that 
A rf will reach 0.9 at T > 0.05t for (n) = 0.94. At the 
temperatures where these calculations were performed, 
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FIG. 5: d-wave irreducible eigenvalue Xd of the Bethe-Salpeter 
equation versus T (a) for U = At and t± = 1.6t, and (b) for 
U = 8t and t± — l.At. Here, the errorbars are smaller than 
the size of the symbols. 

Xd decreases upon doping to (n) = 0.875. Because of the 
QMC "fermion sign problem" [19j, these data were ob- 
tained by using parallel computers. Also shown in Fig. 
5(b) is X d for U = 8t and tj_ = l.At. We find that, at 
T > 0.25i, Xd takes larger values for U = 8t than for 
U = At. These results indicate that the d x 2_ y 2 pairing 
correlations are strong in the Hubbard ladder. 

An important feature of the QMC data presented here 
is that, at sufficiently low T, T s near the Fermi level be- 
comes strongly attractive for q ss momentum transfers. 
This is due to resonant scattering in the d x 2 _ y 2-wa.ve BCS 
channel. In order to demonstrate this effect, consider the 
case of an irreducible interaction T\ which is independent 
of frequency and separable in momentum, 

Wlp) = XV„3 a (p')$«(p) (3) 

OL 

where ev denotes the various pairing channels. In this 
case, the reducible interaction is given by 

r (P'|P) = E 1 y p fl«(p')ga(p), (4) 

a 

with P a = -(T/N)j: p gl(p)\G( P )\ 2 . In general, Tj 
has both attractive and repulsive V a . In Eq. J3J, it is 
seen that the repulsive components get suppressed by re- 
peated scatterings in the BCS channel, while the attrac- 
tive components get enhanced. For the d x 2_ y 2 channel, 
9d(p) — h(cosp x — cospy). If the d x 2_ y 2 component of 
Ti is attractive and gets sufficiently enhanced, here it is 
observed that T for q = p' — p « can become attrac- 
tive. When a d x 2_ y 2-w&ve BCS instability is approached, 
then, at the Fermi surface, r s diverges to — oo for q ss 
scatterings and to +oo for q « (it, it) scatterings. The 
important point is that the exact QMC data for T s also 
exhibit this type of momentum dependence in the Hub- 
bard ladder at temperatures which are not low. 

In summary, we have presented QMC data on the re- 
ducible particle-particle interaction T s in order to de- 



termine the strength and the temperature scale of the 
d x 2_ y 2 pairing correlations in the two- leg Hubbard lad- 
der. We found that T s displays diverging behavior as the 
temperature decreases. In particular, T s exhibits reso- 
nant scattering in the d x 2_ y 2 pairing channel as it would 
be expected for a system near a d x 2_ y 2 BCS instabil- 
ity. We have also shown results on the solution of the 
Bethe-Salpeter equation to quantitatively determine the 
strength of the d x 2_ y 2 pairing. The QMC data presented 
in this paper indicate that the d x 2_ y 2 pairing correlations 
in the Hubbard ladder are strong for certain values of 
the model parameters. We note that, for comparison, 
it would be useful to investigate the strength of d x 2_ y 2 
pairing for a two-leg ladder system which includes both 
onsite Coulomb and electron-phonon interactions. 
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